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We study the precision with which the sign of 5m 2 2 can be determined at a neutrino 
■ factory, as a function of stored muon energies and baselines. This is done by simul- 

taneously fitting the channels Vp — * , v e — > from /i _ decays and the channels 
V[i — > , v e — ► Vp from decays. For a 20 GeV muon storage ring we investigate as 
O ■ a function of the baseline length and the magnitude of 5m\ 2 > the minimum value of the 

parameter sin 2 2#i3 for which the sign of 5m^ 2 can be determined to at least 3 standard 
deviations. We find that for baselines longer than ~ 3000 km, the sign of <5m 2 2 can be 
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determined for sin 26»i3 down to the 10" 3 level with 10 20 decays of 20 GeV muons. 



1 Introduction 



The conceptual development of very intense muon sources has led to a proposal f2| to use 
this new accelerator technology to build a Neutrino Factory in which an intense low energy 
muon beam is rapidly accelerated and injected into a storage ring with long straight sections. 
The muons decaying in these straight sections produce intense beams of highly collimated 
neutrinos. The technical and physics possibilities for building and using neutrino factories 
have been explored by many groups |3|, |J. The interest in neutrino factories is primarily 
driven by recent measurements from the SuperKamiokande (SuperK) collaboration ||, which 
indicate that muon neutrinos produced in atmospheric interactions of cosmic rays oscillate 
into other neutrino flavors, a result that is consistent with measurements made in other 
experiments ||. The presence of a well-defined electron neutrino flux in the muon storage 
ring neutrino beams permits one to explore the effect of matter on the propagation of electron 
neutrinos. The analysis of the oscillation data leads to several scenarios of neutrino masses 
and mixing. In this paper, we investigate the effect of passage through matter on neutrino 
oscillations. We specifically consider the Large Angle MSW scenario (LMA) J7|, which, as 
defined in the recent Fermilab six month physics study ||, has 
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and the CP violating phase 5 = 0. However, to the extent that the contributions from the 
subleading solar Sm^i scale are small, our results apply approximately to other solar sce- 
narios. The transition probabilities in the leading oscillation approximation for propagation 
through matter of constant density are || [H| [11], [12], [13 
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P(u e ^u T ) = c 2 3 sin 2 20™ sin 2 A™, (2) 
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The oscillation arguments are given by 
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where S is given by 
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The amplitude A for v e e forward scattering in matter is given by 



A = 2V2G F N e E u = 1.52 x l(T 4 eV 2 r e p( g/cm 3 )£( GeV) . (7) 

Here Y e is the electron fraction and p(x) is the matter density. For neutrino trajectories 
that pass through the earth's crust, the average density is typically of order 3 gm/cm 3 and 
Y e ~ 0.5. The oscillation probability P{y e — ► u^) is directly proportional to sin 2 26^, which 
is approximately proportional to sin 2 2^ 13 . There is a resonant enhancement for 
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For electron neutrinos, A is positive and the resonance enhancement occurs for positive values 
of 5m| 2 for cos 2#i3 > 0. The reverse is true for electron anti-neutrinos and the enhancement 
occurs for negative values of 5mJ 2 . Thus for a neutrino factory operating with positive 
stored muons (producing a v e beam) one expects an enhanced production of opposite sign 
charged-current events as a result of the oscillation v e — > if 5m\ 2 is positive and vice 



versa for stored negative beams [0, O, O, MJ, 15, 16 



This enhancement is evident in Fig. [TJ, which shows the ratio of v e — > events from \x~ 
decays to v e — > events from /x + decays for 20 GeV muons and a 50 kt detector, assuming 
the oscillation parameters of Eq. ([I]). The results for two other values of 5m\ 2 are also 
presented in Fig. |1| This figure shows that for larger L the ratio of wrong-sign muon events 
is sensitive to the sign of 5m\ 2 . 

The magnitude of 8m\ 2 is determined from the disappearance of muon neutrinos due to 
the oscillation — *> v T , since it can be shown that for the baselines under consideration 
here, the matter effects are small in this channel (see Fig. 2, curves 2 and 3 of Ref. ||). The 
oscillation probabilityP(z/ At — > u T ) can thus be approximated by cos 4 #i 3 sin 2 2 9 23 sin 2 (A ) as 
though it were in vacuum for baselines as far as 4000 km. (See the section "Method" for a 
correction to this approximation for matter effects for longer baselines.) 

2 Method 

In order to extract both the sign and magnitude of brn\ 2 , we simultaneously fit four channels 
(i) — > (ii) v e — > (iii) — > (iv) u e — > v^. The first two channels are measured 
when N /i - decays occur in the storage ring and the second two channels are observed when 
N n + decays occur [|l~7|] . 

We study the extraction of oscillation parameters in two scenarios: 

• A muon storage ring with iV = 10 20 muon decays and a 50 kiloton detector. 

• A muon storage ring with iV = 10 19 muon decays and a 50 kiloton detector. A 20 GeV 
version of this is known as an "entry level" neutrino factory. 

Our study is performed for baselines ranging from 732 km up to 7332 km and stored 
muon energies ranging from 20 GeV to 50 GeV. We calculate the neutrino event rates by 
propagating the neutrinos through matter, taking into account the variations in the density 



profile using the Preliminary Reference Earth Model [18| by solving the evolution equations 
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numerically. Realistic detector resolutions are used that are appropriate for a magnetized 
iron scintillator detector ]12|, ||. We assume a muon energy resolution J- = 0.05 and a 



hadronic shower resolution of J- = 0.53/ y/Eh for showers of energy Eh > 3 GeV and 
0.8/y/Eh for showers of energy E h < 3 GeV. We accept only events that possess muons of 
true energy greater than 4 GeV. Figure ^| shows the wrong-sign muon appearance spectra 
with these cuts as function of 5m\ 2 for both fi + and fi~ beams for both signs of 5m\ 2 at a 
baseline of 2800 km. The resonance enhancement in wrong sign muon production is clearly 
seen in Fig. [| (b) and (c). Using these histograms and similar ones for the disappearance 
channels, it is possible to predict the spectrum in any channel for any value of 5m 2 2 in the 
range of interest by polynomial interpolation of the histograms bin by bin with the method of 
divided differences. For the interpolation, the disappearance probabilities can be treated to 
first order as though they are due to vacuum oscillations (i.e. matter effects can be neglected) 
and are proportional to cos 4 #13 sin 2 2^23 . For the longer baselines, the matter effects affect 
these disappearance rates slightly, as can be seen from Fig. 3 of Ref . . The difference in the 
disappearance rates for positive and negative 5m\ 2 is proportional to sin 2 20i3, vanishing as 
sin 2 2^i3 — > 0. For the channel — > z/ M , the positive 5m\ 2 solution approaches the negative 
5m| 2 solution as sin 2 26^ — > 0, the reverse being the case for anti- neutrinos. We use this 
behavior to interpolate the disappearance spectra as a function of sin 2 2^13. The event 
rates in both the appearance and disappearance channels depend only on the three vacuum 
parameters sin 2 2^23, sin 2 26 l i3, and 5m| 2 . Using this interpolation scheme, we can generate 
events for any combination of these parameters for various muon momenta and baselines. 
Most of the information on the sign of <5m 2 2 comes from the appearance channels and the 



precision on the magnitude of 5m^ 2 comes from the disappearance channels || 112, [13], [15 



2.1 Backgrounds 

The background for the wrong-sign muon appearance channels arise from three sources: 
(i) charm production, (ii) 7r, K decay events producing wrong sign muons in neutral current 
interactions, and (iii) 71", K decay events producing wrong sign muons in charged-current (CC) 



events where the primary muon was considered lost || |12| . A more detailed study of the 
backgrounds has shown that significant reductions are obtained by demanding that P 2 > 
2 GeV 2 , where P 2 is defined as the transverse-momentum-squared of the muon with respect 
to the hadronic shower direction. Imposing this cut results in a further signal efficiency 
factor of 0.9 |19| for the CC disappearance channels, 0.62 for the /i + beam appearance signal 



and 0.45 for the \i~ beam appearance signal. This results in an average background/V^ — > 
CC signal rate of 0.45 x 10~ 4 for fi + beam appearance channel and 0.25 x 10~ 5 for the \x~ 
beam appearance channel. The difference in these rates between the fi + and \x~ beams is due 
to the different kinematical distributions of the neutrino and anti-neutrino CC interactions. 
We fold these rates into the theoretical prediction for each channel. We also assume a 
normalization systematic uncertainty of 1% between the fi + beam and \i~ beam events. 

Table [1] lists the wrong sign muon appearance- and background-rates for an "entry level" 
machine. Event rates for the appearance channels are appreciable, whereas background rates 
are negligible for baselines longer than about 2800 km. 
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2.2 Fitting 



Events are generated for the 4 channels — > u^, v e — > z/ M , ^ — ► z/ M , and i> e — > and for 
both signs of 5m| 2 - The four simulated measured event rates are fitted simultaneously for the 
three parameters sin 2 2#i3, sin 2 2^23, and 5m 2 2 . In each fit the sign of 5m\ 2 is constrained [pO 



The difference in negative log-likelihood, AL, between the fits in which sign of 5m| 2 has been 
constrained to the correct and incorrect sign is evaluated for a given value of sin 2 2# 13 . This 
is then repeated for a range of sin 2 26^. It is empirically found that the average AL varies 
linearly with sin 2 2^13 and this is used to estimate that value of sin 2 2^13 at which AL=4.5, 
i.e a Gaussian 3a ability to differentiate the sign of 8m\ 2 . This value of sin 2 2^13 we define 
as the 3a reach in sin 2 2#i3 space. We only fit for the number of events and not the shape 
of the spectra, since at the 3a point statistics are such that shape information contributes 
little. 



Table 1: Wrong-sign muon rates for a 50 kt detector (with a muon threshold of 4 GeV) a 
distance L downstream of a muon factory (energy E^) providing 10 19 muon decays. Rates 
are shown for LMA scenario of Eq. (j^) with both signs of 8m\ 2 considered separately. The 
background rates listed are for each sign of 5m\ 2 and do not depend on the sign of <5m 2 2 . 
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3 Results 

Figure [3] shows the difference in negative log-likelihood between a correct and wrong-sign 
mass hypothesis expressed as a number of equivalent Gaussian standard deviations versus 
baseline length for muon storage ring energies of 20, 30, 40 and 50 GeV. The values of the 
oscillation parameters are for the LMA scenario in Eq. [I]. Figure 0(a) is for 10 20 decays 
for each sign of stored energy and a 50 kiloton detector and positive 5m\ 2 (b) for negative 
5m\ 2 for various values of stored muon energy. Figures |3| (c) and (d) show the corresponding 
curves for 10 19 decays and a 50 kiloton detector. An entry-level machine would permit one 
to perform a 5a differentiation of the sign of 5m| 2 at a baseline length of ~2800 km. 
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Figure [| shows the 3cr reach of sin 2 2# 13 versus baseline length for a 20 GeV muon storage 
ring (a) with 10 20 decays, a 50 kiloton detector and positive 5m\ 2l (b) for negative 5m\ 2 for 
various values of |5mg 2 |. Figures |] (c) and (d) show the corresponding curves for 10 19 
decays and a 50 kiloton detector. The error bars show the uncertainties due to statistical 
fluctuations in our determination of the 3a reach. Our results agree with similar calculations 



in Ref. [15]. It can be seen that an entry-level machine is capable of determining the sign of 



5r77,3 2 provided that sin 29 13 is greater than 0.01 for bm\ 2 m the range 0.0025-0.0045 eV . 

To conclude, we have shown that the neutrino factory provides a powerful tool to explore 
the structure of neutrino masses. Even an "entry level" machine with 10 19 decays of 20 
GeV stored muons and a 50 kiloton detector at a baseline exceeding ~2800 km is capable of 
differentiating the sign of 8m\ 2 for the LMA scenario. The 3cr reach in sin 2 26^ in the sign 
determination of 8m\ 2 for such a machine exceeds the present bound on sin 2 2# 13 by about 
an order of magnitude. 
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Figure 1: The ratio of wrong-sign muon event rates N(u e — > v^)jN{y e — > v^) versus baseline 
for a 20 GeV muon storage ring and 8m\ 2 = 0.0025 eV 2 (dashed line), 0.0035 eV 2 (solid line), 
and 0.0045 eV 2 (dotted line). The other oscillation parameters are given in Eq. (|l|). A 4 GeV 
minimum cut was imposed on the detected muon energy. The error bars show the statistical 
errors corresponding to 10 20 decays and a 50 kiloton detector. 
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Figure 2: The wrong sign rnuon appearance rates for a 20 GeV rnuon storage ring at 
a baseline of 2800 km with 10 20 decays and a 50 kiloton detector for (a) fi + stored and 
negative 5m\ 2 , (b) /i~ stored and negative 5m\ 2 , (c) fi + stored and positive 5m\ 2 , (d) 
stored and positive 8m\ 2 - The values of |<5t«3 2 | range from 0.0005 to 0.0050 eV 2 in steps of 
0.0005 eV 2 . Matter enhancements are evident in (b) and (c). 
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Figure 3: The statistical significance (number of standard deviations) with which the sign 
of 5m^ 2 can be determined versus baseline length for various muon storage ring energies. 
The results are shown for a 50 kiloton detector, and (a) 10 20 fi + and fi~ decays and positive 
values of 8m\ 2 ', (b) 10 20 ji + and pT decays and negative values of <5m 2 2 ; (c) 10 19 /i + and \i~ 
decays and positive values of 8m\ 2 ] (d) 10 19 /i + and pT decays and negative values of 8m\ 2 . 
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Figure 4: The minimum value of sin 2 26i% for which the sign of 5m\ 2 can be determined to 
at least 3 standard deviations versus baseline length for various values of 5m\ 2 . The results 
are shown for a 50 kiloton detector, and (a) 10 20 /i + and fi~ decays and positive values of 
Sm 2 32 ; 



(b) 10 20 fi + and \i decays and negative values of 5m\ 2 ] (c) 10 19 /i + and ji decays 



and positive values of 5m| 2 ; (d) 10 19 /i + and \i decays and negative values of 8m\ 2 . 
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